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Abstract. We present the results of an investigation of post-Main Sequence mass loss from stars in clusters in the 
Magellanic Clouds, based around an imaging survey in the L'-band (3.8 pm) performed with the VLT at ESO. 
The data are complemented with JHK S (ESO and 2MASS) and mid-IR photometry (TIMMI2 at ESO, ISOCAM 
on-board ISO, and data from IRAS and MSX). The goal is to determine the influence of initial metallicity and 
initial mass on the mass loss and evolution during the latest stages of stellar evolution. Dust-enshrouded giants 
are identified by their reddened near-IR colours and thermal-IR dust excess emission. Most of these objects are 
Asymptotic Giant Branch (AGB) carbon stars in intermediate-age clusters, with progenitor masses between 1.3 
and ~5 Mq. Red supergiants with circumstellar dust envelopes are found in young clusters, and have progenitor 
masses between 13 and 20 Mq. Post- AGB objects (e.g., Planetary Nebulae) and massive stars with detached 
envelopes and/or hot central stars are found in several clusters. We model the spectral energy distributions of 
the cluster IR objects, in order to estimate their bolometric luminosities and mass-loss rates. The IR objects are 
the most luminous cluster objects, and have luminosities as expected for their initial mass and metallicity. They 
experience mass-loss rates in the range from a few 10 -6 up to 10 -4 Mq yr -1 (or more), with most of the spread 
being due to evolutionary effects and only a weak dependence on progenitor mass and/or initial metallicity. About 
half of the mass lost by 1.3-3 Mq stars is shed during the superwind phase, which lasts of order 10 5 yr. Objects 
with detached shells are found to have experienced the highest mass-loss rates, and are therefore interpreted as 
post-superwind objects. We also propose a simple method to measure the cluster mass from L'-band images. 

Key words. Stars: AGB and post-AGB - Stars: evolution - Stars: mass-loss - supergiants - Magellanic Clouds - 
Infrared: stars 
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carbon stars are more luminous than the M stars, with the 
S stars in between. The transition luminosity M-S and S- 
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stellar dust emission are essential to measure the mass-loss 
rate from these stars. 

Only a few cluster IR objects are known. These in- 
clude the Mi ~ 4 Mq OH/IR star IRAS 05298—6957 
(van Loon et al. 2001a) and the M; ~ 2.2 M Q carbon 
star LI-LMC1813 (van Loon ct al. 2003). Tanabe et al. 
(2004) surveyed clusters in the Magellanic Clouds with 
the ISOCAM instrument onboard ISO. They detected 
two dust-enshrouded AGB stars in each of the populous 
intermediate-age clusters NGC419 (SMC) and NGC 1783 
and NGC 1978 (LMC), four of which were already known 
from near-IR data. 

We here present results from a systematic survey for 
dust-enshrouded red giants in clusters in the Magellanic 
Clouds. The investigation is built around L'-band imag- 
ing obtained at the ESO-VLT, supplemented by near-IR 
imaging photometry and a collection of mid-IR data from 
IRAS, ISO, MSX and the ESO 3.6m telescope. We find 
30 dusty stars, of which 26 are likely cluster members. 
We determine their nature, luminosity and mass-loss rate, 
and investigate the dependence on progenitor mass and 
metallicity. 
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Fig. 1. Age- metallicity diagram for the SMC (open sym- 
bols) and LMC (solid symbols) clusters in our sample. 



2. Cluster selection, properties, and bias 

We have attempted to select clusters across a wide range 
of ages and metallicities (Table 1, criterion "C"). Although 
the more populous clusters are attractive targets to sur- 
vey for rare objects such as dust-enshrouded stars, small 
clusters of t^l Gyr are much more numerous and many of 
the IR-object selected clusters (Table 1, criterion "I") turn 
out to be associated with this type of cluster. The draw- 
back of these less conspicuous clusters is that they have 
not usually been studied in detail and ages and metal- 
licities are often inaccurate — if at all available. Due to 
the smaller stellar over-density with respect to the field, 
membership of small clusters is also more ambiguous. We 
added from the literature (Table 1, criterion "L") another 
two clusters in which IR objects had been found. 

The cluster ages and metallicities — for which we take 
the commonly used relative iron abundance [Fe/H] — are 
collected from a variety of recent sources in the literature, 
and are listed in Table 1. For some clusters, average values 
were taken from multiple sources. In the case of NGC 1978 
the disagreement between the individual sources is rather 
large, with Olszewski et al. (1991) and Hill et al. (2000) 
disagreeing over the metallicity by 0.5 dex for the spec- 
troscopic determinations of two stars in common. 

The clusters show a clear correlation between age 
and metallicity (Fig. 1), which is observed for both the 
SMC and the LMC (e.g., Olszewski et al. 1991; Geisler 
et al. 1997; Da Costa & Hatzidimitriou 1998; de Freitas 
Pacheco, Barbuy & Idiart 1998). The turnover to lower 
metallicities of the youngest clusters is intriguing, with 
young, fairly massive clusters as metal-poor as [Fe/H] = 
-0.9 observed in both clouds (Hill ct al. 2000). There is 
sufficient spread in the cluster properties that one cannot 



assume that every LMC cluster is more metal-rich than 
an SMC cluster of the same age. 

Mackey & Gilmore (2003a, b) estimated the cluster 
masses for 10 clusters in the SMC and 53 clusters in the 
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Fig. 2. Mass-radius diagram for clusters in the SMC 
(open symbols) and LMC (solid symbols), with masses 
from Mackey & Gilmore (2003a,b) and radii from Bica & 
Schmitt (1995) and Bica et al. (1999). The bottom panel 
shows the mass ratio with respect to the displayed param- 
eterization M r t . 
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Table 1. List of programme clusters, in order of increasing Right Ascension (all coordinates arc in J2000). The cluster 
radii are taken from Bica & Schmitt (1995) and Bica et al. (1999) for the SMC and LMC, respectively, and converted to 
parsecs assuming distances to the SMC and LMC of 60 and 50 kpc, respectively. Metallicities without error estimates 
are assumed values; masses without error estimates are derived from the mass-radius-age relation as described in 
the text. Also given are the numbers of stars detected in the L'-band images with an L'-band brightness exceeding 
that of the RGB tip, at Mu < —6.4 mag (see Section 5.5). The references for the ages, metallicities and masses are 
as follows: l=Beasley, Hoyle & Sharpies (2002); 2=Bica et al. (1996); 3=de Freitas Pacheco et al. (1998); 4=Elson 
(1991); 5=Kontizas (1980); 6=Mackey & Gilmore (2003a); 7=Mackey & Gilmore (2003b); 8=Mateo (1988); 9=Mighell, 
Sarajedini & French (1998); 10=Oliva & Origlia (1998); ll^Olszewski et al. (1991); 12=Pictrzyriski & Udalski (2000); 
13=Vallenari, Bettoni & Chiosi (1998); 14=van Loon et al. (2001a); 15=van Loon et al. (2003); x=average of de Freitas 
Pacheco et al. (1998), Hill et al. (2000) & Olszewski et al. (1991); y=this work. The last column indicates the criterion 
upon which the selection was based (C=" cluster" , I="IR", L= "Literature"). 
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LMC, some of which are in common with our sample, weak correlation between the mass and volume (Fig. 2, 
We collected the radii for all their clusters from Bica & top panel): 
Schmitt (1995) and Bica et al. (1999). We thus find a 



log M r [M ] = 1.12 + 31ogr[pc], 



(1) 
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with a standard deviation of 0.7 in logM. The scatter 
can be reduced by a parameterization in both radius and 
cluster age (Fig. 2, bottom panel): 

logM r , t [M ] =a + &logr[pc] + clogi[yr], (2) 

where a = -0.7±0.5, b = 1.3±0.4 and c = 0.44±0.06. We 
use the M r j relation to estimate the mass of clusters that 
were not included in the analysis of Mackey & Gilmore 
(Table 1). 

Throughout this paper we use distance moduli to the 
SMC and LMC of 19.0 and 18.5 mag, respectively. 

3. Observations 

3.1. J and K s -band imaging 

Images were obtained through the J (1.25 fim) and K s (2.2 
/im) bands with the SOFI instrument on the ESO 3.5m 
NTT at La Silla, Chile, on the first half of the night of 7 
November 2000. At a scale of 0.146" pixer 1 , the field-of- 
view was 2.5' x 2.5'. The images were constructed through 
combination of 9 frames which were obtained at slightly 
different (few arcsec) "jitter" positions to allow for the 
construction and subtraction of the (bright and variable) 
sky background. The seeing was typically around 1", but 
some of the images were distorted because the delay time 
to allow the telescope to recover from its movement to a 
new jitter position was set too short. The weather condi- 
tions deteriorated shortly after completion of the cluster- 
selected subsample, and consequently the J and K s -band 
data obtained for the IR-object-selected subsample are of- 
ten quite poor and in several cases not obtained at all. The 
total integration time per final image is 90 seconds except 
where images were discarded because no stars were visible 
due to cloud cover. 

The data were reduced using the ESO Eclipse and 
MIDAS software packages. First, electronic ghosts were 
removed with help of the Eclipse tool IS.GHOST. Then, 
a flatfield frame was constructed from the differences be- 
tween the pair-wise observations of an illuminated and 
dark screen, to correct for transmission, vignetting and 
detector response variations across the frame. Next, each 
frame was sky-subtracted by first (for each pixel) subtract- 
ing the minimum value of the set of 9 jittered frames (after 
having scaled each frame to its median value), and then 
subtracting the pixel value obtained by median-averaging 
along the columns of the frame. The resulting frames were 
shift-added to create the final image. 

Photometry was obtained for all the stars detected 
in the L'-band images described below, using (software) 
apertures with a diameter of 1.6". A few stars in NGC 419, 
NGC1984 and NGC 2210 were measured through 0.7"- 
diameter apertures in the case of exceptional crowding, 
and their magnitudes were corrected for the aperture dif- 
ference. The magnitudes were calibrated against those of 
known 2MASS sources in the images to a systematic ac- 
curacy of ~0.1 mag, but relative photometry within the 
images themselves is generally better. 



2MASS data are used where ESO J and K s -band data 
are absent, as well as to complement the analysis with 
H-band (1.65 /im) photometry. The J-band magnitude of 
NGC 1994-IR1 and the K-band magnitudes of two other 
bright stars in the same cluster were reconstructed by dif- 
ferential aperture photometry on the 2MASS image. The 
2MASS data are problematic for clusters in the Magellanic 
Clouds though, as many 2MASS sources are in fact blends 
of stars. Also, the 2MASS data were obtained more than 
two years before the L'-band data whilst the ESO J and 
K s -band data were obtained within a few weeks after the 
L'-band data. Many of the dust-enshrouded red giants are 
large-amplitude variables with periods of around 400-1000 
days and hence, where the ESO photometry can be consid- 
ered near-simultaneous with the L'-band photometry, the 
2MASS photometry have random phase differences with 
respect to the L'-band photometry. 

Some objects have been detected in the DENIS survey 
of the Magellanic Clouds (Cioni ct al. 2000). Although 
less deep than 2MASS, it provides an I-band magnitude 
for some of the (optically) brighter, relatively isolated ob- 
jects. This is useful in constraining the spectral energy 
distributions of hot stars such as post-AGB objects and 
PNe. The DENIS catalogue also lists B and/or R-band 
magnitudes in some cases. Crowding is more severe in the 
optical, though, which could render the optical colours of 
intrinsically red objects bluer due to the contribution at 
short wavelengths of other, less red stars. 

3.2. L'-band imaging 

Images were obtained through the L' band (3.78 ^m, 
bandwidth AA = 0.58 ^m) with the ISAAC instrument on 
the ESO VLT at Paranal, Chile, in service mode during a 
total of three hours spread over the nights of 12, 17 and 
18 October 2000. At a scale of 0.071" pixel" 1 , the field-of- 
view was 72" x 72". The images were obtained following 
the usual chop-and-nod procedure for thermal-IR obser- 
vations, where both the chop and the nod throw were 10" 
in the North-South direction. The total exposure time per 
target is two minutes, with a detector integration time of 
0.104 sec. The conditions were photometric, with a seeing 
of typically around 0.5". 

The data were reduced using the ESO Eclipse and 
MIDAS software packages. First, electronic ghosts were re- 
moved with help of the Eclipse tool IS.GHOST. A flatfield 
frame was constructed by imaging the twilight sky and 
using the Eclipse tool is_twflat. After flatfielding the 
science frames obtained at each nod position, the result- 
ing pair of frames were subtracted and then shift-added to 
produce a frame with stellar images that contain all of the 
observed light. Photometric measurements were obtained 
through a 1.07" software aperture, and were calibrated 
against the measurements of the standard stars HD 22686, 
HD 19904 and HR2354. 

Additional photometry was obtained from acquisition 
images for L-band spectroscopy for the clusters NGC 1903, 
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SL 453, NGC 1984, NGC 1978 and SL 519 (7 & 8 December 
2003; van Loon et al., in preparation) and for KMHK 1603 
(14 December 2001; van Loon et al. 2003). These im- 
ages were obtained through the L'-band filter except for 
NGC 1978 and the standard star HIP 020020 A+B that 
were observed through narrow-band filters centred at 3.8 
and 3.28 /zm, respectively. Detector integration times were 
0.025 sec for the clusters and 0.28 sec for NGC 1978 and 
the standard star. The same zeropoints were assumed for 
all images, which by cross-correlation with our dedicated 
L'-band photometry seemed to be accurate to within a 
few 0.1 mag. 

3.3. N-band imaging 

Images were obtained in the N band (8-13 /im) with 
the TIMMI2 instrument on the ESO 3.6m at La Silla, 
Chile, on the nights of 19-21 January 2001, of the clusters 
NGC 330 and KMHK 1603 through the N2 filter (central 
wavelength 10.6 /im, 9.7-11.3 /im at 50% transmission) 
and of the clusters NGC 419, NGC 1984 and NGC 1994 
through the 11.9 /im filter (11.0-12.2 /im at 50% trans- 
mission). At a scale of 0.3" pixel -1 (0.2" pixel -1 for 
NGC 330), the field-of-view was 96" x 72" (64" x 48" 
for NGC 330). Standard chopping and nodding procedures 
were followed, with a throw of 10". Total effective expo- 
sure times amounted to 4-9 minutes per cluster, using 
detector integration times of 22 to 27 msec. 

The data were reduced using the ESO MIDAS soft- 
ware package. After inspection of all individual frames 
and rejection of saturated data, the frames belonging to 
the same nod position were averaged and then subtracted 
from the average frame at the complementary nod posi- 
tion. The thus obtained images were shift-added to con- 
struct the final image, from which by additional shift- 
adding an image was produced with stellar images that 
contain all the observed light. On these stellar images 
photometry was obtained through a 3" (software) aper- 
ture and which was calibrated against measurements of 
the standard stars A2 Tuc, aHya, 7 Pic, 7 Ret and 9 Dor. 
An early result for LI-LMC 1813 (KMHK 1603-IR1) has 
already been presented in van Loon et al. (2003). 

3.4. Additional mid-1 R imaging 

For the identified cluster IR objects, additional mid-IR 
data was obtained from the IRAS and MSX surveys and 
in some cases from ISO survey data. Like 2MASS, these 
data are not contemporary with the L'-band data, how- 
ever they are important in constraining the spectral en- 
ergy distribution to estimate the bolometric luminosity 
and mass-loss rate. 

For each IR object we collected the corresponding 
IRAS scans from the IRAS data server 1 and used the 
GIPSY software with the SCANAID tool to reconstruct a 
cut through the emission on the exact position. This is pre- 

1 http:/ /www. astro. rug.nl/IRAS-Server/ 



ferred over using the IRAS Point Source Catalogue values 
and works well for isolated sources, for which often reliable 
flux densities can be estimated down to a level of a few 
0.01 Jy at 12 and 25 /zm and ^0.1 Jy at 60 /tm, but where 
clusters contain more than one IR-bright source the flux 
densities can become quite unreliable and should often be 
regarded as upper limits. 

Most cluster IR objects have reliable MSX band A (8.3 
/im) flux densities in the MSX Point Source Catalogue 
(Version 2.3). However, the flux densities for NGC 419- 
IR1 and IR2, NGC 1903-IR1 and IR2, NGC 1978-IR3 and 
IR4, HS327-E-IR2, SL519-IR2 and NGC2121-IR1 were 
estimated by differential aperture photometry on the MSX 
images with respect to nearby moderately-bright sources 
listed in the MSX Point Source Catalogue. The SMC clus- 
ter objects NGC 419- IR1 and IR2 were also measured on 
the original MSX images. Only a few objects are bright 
enough to have reliable MSX measurements in bands C 
(12.1 /im), D (14.7 itm) or E (21.3 /zm). 

Three of the clusters, BSDL 1837, HS327 and SL519 
were covered by the ISO mini-survey of the Magellanic 
Clouds (Loup et al., in preparation), performed with the 
ISOCAM instrument and filters LW1 (4.5 /im) and LW10 
(12.0 fim). SL519 is located in the prime calibration field 
for the survey, which is also one of the calibration fields 
for the ISOGAL survey (Omont et al. 2003), and was 
therefore also observed twice with the LW2 (6.7 /Km) fil- 
ter, a second time with the LW10 filter and once with 
the LW3 (14.3 /im) filter. Flux densities were measured 
through differential aperture photometry with respect to 
moderately-bright objects in the ISO mini-survey and/or 
ISOGAL Point Source Catalogues. Time-averaged images 
were created for SL519 by rotate-shift-adding the images 
taken at the different epochs. 

4. Results 

4.1. Identification of IR objects 

The cluster IR objects and their identification with known 
mid-IR sources, IR photometry and classification, and the 
(expected) availability of 3-4 /zm spectroscopy and/or 
Spitzer Space Telescope IRS observations are listed in 
Tables 2-7. Table 2 includes literature references for the 
individual IR objects, on a few of which we comment here: 

- The variable Ml. 5 supergiant HV 12501 (WOHS78) 
has been held responsible for the mid-IR emission in 
KMHK 292 (Trams ct al. 1999b). However, another 
supergiant in the cluster, WOHS76 is brighter than 
HV 12501 both in the ncar-IR and in MSX band A. 

- NGC 1852-IR1 is identified with the low-excitation PN 
SMPLMC31 (Morgan 1984). 

- NGC 1984 is associated with the bright IR source 
IRAS 05280-6910 as well as a source of OH (Wood 
et al. 1992) and H 2 (van Loon et al. 2001b) maser 
emission. ATCA observations (mentioned by Wood et 
al. 1992, but unpublished) locate the OH maser with 
RA= 5 h 27 m 39.87 s and Dcc= -69°08'06.8" (J2000) 
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Table 2. List of cluster IR objects (J2000 coordinates, based on 2MASS), and identifications with mid-IR sources. 
Mid-IR flux densities (in Jy) were obtained from the MSX Point Source Catalogue or from the original MSX images, 
and from IRAS scans; the bands are listed by their central wavelengths in fxm. Values marked with a colon are suspect. 
The references are as follows: l=Egan, Van Dyk & Price (2001); 2=Elias, Frogcl & Humphreys (1985); 3=Ferraro et 
al. (1995); 4=Frogel et al. (1990); 5=Israel & Koornnecf (1991); 6=Leisy et al. (1997); 7=Lloyd-Evans (1980); 8=Loup 
et al. (1997); 9=Marshall et al. (2004); 10=Matsuura et al. (2002); ll=Morgan (1984); 12=Nishida et al. (2000); 
13=Reid, Tinney & Mould (1990); 14=Sanduleak, MacConnell & Philip (1978); 15=Tanabe et al. (1997); 16=Tanabe 
et al. (1999); 17=Tanabe et al. (2004); 18=Trams et al. (1999a); 19=Trams et al. (1999b); 20=van Loon, Zijlstra & 
Groenewegen (1999a); 21=van Loon et al. (1997); 22=van Loon et al. (1998); 23=van Loon et al. (1999b); 24=van 
Loon et al. (2001a); 25=van Loon et al. (2001b); 26=van Loon et al. (2003); 27=van Loon et al. (2005); 28=van Loon et 
al., in preparation; 29=Villaver, Stanghellini & Shaw (2003); 30= Westerlund, Olander & Hedin (1981); 31=Whitclock 
et al. (2003); 32=Wood & Cohen (2001); 33=Wood et al. (1992). 



Name 






RA 


L/cL 


Identification(s) 


1V1D-A. 


TDAC 
irv.iT-u 








Cluster IR# 




m 


\ ) 


TD AC 

1XX-TVO 


IVloTv 


o.O 


1 9 
1Z 


9^ 
ZO 


fin 
ou 




bmall Magelt 


anic 


Cloud 


















NGC 419 


1 


1 


08 12.97 


79 E.9 A A fl 
— ( Z Oi 44. U 






U.UOo 


<U.1U 


^(\ 1 A 

<U.14 




19 1 ^ 17 

lz,10,l r 




2 


1 08 17.52 


79 E.Q no O 

-iz oo uy.z 


-Lt-olVlO loz 




U.loo 


n oa 
U.Z4 


n in 

u.iy 


<u.z 


if; i 7 
10,1 / 


Large Magell 


anic 


Cloud 


















TJG OO 

rib 33 


1 


4 


AC\ 1 O A f\ 

49 18. 4y 


— oy oo 14.0 


U44yo — OyDo 


lloU 


n A C\A 
U.4U4 


U.ol 


n 99 a 
U.zz 


<U.l 


i o in on oo oq oi 

lo,ly,zU,zz,zo,ol 


rvlvlrlrv 29z 


1 


4 


55 34.85 


fiG 9fi E.E. 1 

— Oy zO 00. / 


U4ooy— Dyol 


1 Q9G 

lozy 


U.zyo 


n Q9 
U.oz 


U.lo 


<U.z 


oU 




o 
z 


4 


ob 41. oo 


gs\ Ofi OA Q 

— oy zo zq.o 






U. 1 ( 1 


n 90- 
u.zu. 


u.uo. 




1 Q 9Q Qfl 

iy,zo,ou 


rvlvlrlrv zoo 


-1 

1 


A 
4 


EE oo no 
00 Oo.yo 


an /in in 7 
— D / 4y 1U. ( 


U400 / — O / Do 


1 9Q8 

izoo 


U. 10 ( 


n oa 

U.Z4 


U. 10 


^U.O 


in 1Q 99 Qi 

iu,iy,zz,oi 


AT*^ f ' 1 TOO 

JNGC 1783 


1 


4 


59 01.11 


gx: co on Q 
— DO Oo OU.O 




1 97*3 
1Z 1 


U. 1 10 


u.uo 


u.uo 




19 1^ Ifi 17 
1Z, 10, 10,1 / 


IN GO looz 


1 


5 


09 20.22 


fi7 Al 9^ n 
— 01 4/ ZO.U 


uouy4 — d / or 


oo 


U.UOO 


U. 10 


U.40 




fill 1Q 1/1 OQ 

0,11, io, 14, zy 


IN GO 19U3 


-1 

1 


5 


17 16.33 


en on on o 

— oy zu zy.o 






U. 1U4 


9n 
<^u.zu 


9n 
<^u.zu 




9ft 
Zo 




o 
z 


5 


17 17.38 


fiG on ^A fi 
— oy zu o4.o 






u.uoz 


90 

<^u.zu 


90 
<^u.zu 




98 
Zo 




3 


5 


17 22.62 


fiQ on 1 ^ ^ 
— oy zu 10.0 




344 


U.OUl 


n 99- 
u.zz. 


n no- 
u.uy. 






HS270 


1 


5 


23 53.93 


—71 ^4 4S Q 


05246-7137 


423 


0.161 


0.19 


0.40 


2.0 




SL453 


1 


5 


25 03.26 


-69 26 17.0 




484 


0.156 


<0.10 


<0.10 




28 


SL482 


1 


5 


27 17.84 


-66 22 05.6 


05273-6624 


562 


0.233 


<0.10 


<0.12 




8,13 


NGC 1984 


1 


5 


27 40.83 


-69 08 05.4 


05280-6910 




£1.370 


3.90 


23.50 


12.6 


9,25,28,32,33 




2 


5 


27 40.11 


-69 08 04.5 






<1.370 








2,28,30 




3 


5 


27 35.67 


-69 08 56.3 






0.157 


0.10: 


<1.00 




5,6,11,14 


BSDL 1837 


1 


5 


27 47.48 


-69 13 20.5 


05281-6915 


588 


0.320 


0.30 


0.20 




30 


NGC 1994 


1 


5 


28 21.98 


-69 08 33.7 


05287-6910 


583 


£0.962 


1.07 


1.40 


<1.0 




NGC 1978 


1 


5 


28 40.17 


-66 13 54.2 






0.077 


0.07: 


<0.05 




4,12,17,28 




2 


5 


28 47.20 


-66 14 13.6 


05287-6616 




0.124 


0.08: 


<0.10 




3,16,17,28 




3 


5 


29 02.41 


-66 15 27.8 


05289-6617 




0.080 


0.16 


0.38 


0.3: 


13,19,23,27 




4 


5 


28 44.50 


-66 14 04.0 






0.070: 


0.10: 


0.10: 




4,7,17 


HS327-E 


1 


5 


29 24.60 


-69 55 13.4 


05298-6957 


653 


0.292 


0.85 b 


1.38 b 


<3.0 


9,19,22,23,24,25,33 




2 


5 


29 25.46 


-69 54 52.3 






0.030: 


<0.85 






24 


SL519 


1 


5 


30 04.87 


-69 56 45.2 




654 


0.152 


0.20: 


0.10: 




28 




2 


5 


30 02.15 


-69 56 17.2 






£0.030 


0.15: 


0.10: 






NGC 2100 


1 


5 


42 11.57 


-69 12 48.8 


05425-6914 


1435 


0.167 


0.30 


0.22 




1 


NGC 2121 


1 


5 


48 16.81 


-71 28 39.3 






0.028: 


0.02: 


<0.05 




7,4,17 


KMHK 1603 


1 


6 


02 31.06 


-67 12 47.0 


06025-6712 


1652 


0.251 


0.39 


0.28 


£0.2 


21,26 
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very near the Ml red giant WOH G 347. However, our 
L'-band images reveal a very red object at only a 
few arcseconds from WOH G 347. Our mid-IR image 
taken with TIMMI2 at a wavelength of A = 11.9 /im 
clearly identifies this red star with the IRAS source: 
WOH G 347 contributes only 5 per cent to the total 
emission at this wavelength. The IRAS source and the 
cluster itself have been mis-identified in the past (e.g, 
in Simbad) with SMP LMC 64, a very low excitation 
PN almost an arcmin away. 



— SL519-IR1 dominates the IR emission seen in the 
ISOCAM images (Fig. 3). Another very faint, some- 
what red object is visible too: SL519-IR2. There are 
three other not very red mid-IR sources which however 
are each blends of at least 2-3 bright individual stars. 

The cumulative near-IR colour-(absolute)magnitude 
diagrams (Fig. 4) of all L'-band detected stars in and 
around all observed clusters in the SMC (open symbols) 
and LMC (dots) are quite similar for both Magellanic 
Clouds. They show the prominent branches of RGB (up 
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Table 3. Ncar-IR photometry of cluster IR objects: JHK S magnitudes from 2MASS and/or JK S L' magnitudes from 
our own imaging ( "ESO" , and "acq" for the spectroscopy acquisition) . 
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to Mks — —6 mag), AGB (up to Mk s — —9 mag), mas- 
sive main-sequence stars and RSGs (the latter dominating 
the IR light, at Mk s Ss — 10 mag). In the J-K s colour, a 
separate branch of massive hot stars stands out, which is 
however difficult to distinguish in the K s -L' colours. At 
first glance it seems that the SMC supergiants are bluer 
than similar LMC supergiants by about 0.2 mag in both 
J-K s and K s -L', yet the upper AGB sequences seem to 
coincide in the SMC and LMC. 

The bright (Mk s < —10 mag) stars with negative 
K s -L' colours are all from the clusters NGC 1994 and 
NGC 2004 for which only 2MASS K s -band magnitudes are 
available: crowding is a severe problem for these compact 
clusters, and several of the brightest cluster members be- 
come blended in the 2MASS data leading to overestimates 
in the individual stars' K s -band magnitudes. The L'-band 
data do not suffer from the crowding and hence the K s -L' 
colours are being under- estimated. Because the blending 
of similar stars does not affect their 2MASS colours the 
effect is not apparent in the J-K s colours. 

Stars with colours (J — K s )^2 or (K s — L')^l mag are 
considered to be "IR objects". These have been reddened 



compared to normal photospheric colours by circumstel- 
lar or interstellar selective extinction (notably affecting 
the J-K s colour) or by circumstellar emission giving rise 
to excess emission at thermal IR wavelengths (e.g., L'). 
It is worth pointing out that the reddest SMC object, 
NGC419-IR1 is as red as the reddest LMC objects. All 
objects with (K s — L') > 3 mag and Mk s > —6 mag are 
undetected in the J-band. This includes four out of five 
objects that form a clump in the colour-magnitude dia- 
gram around (K s — L') ~ 3.4, Mk s — —5.6 mag. A few 
objects are included in our analysis that are not very red 
but dominate the cluster's light at mid-IR wavelengths. 

4.2. IR colour-magnitude diagrams of selected clusters 

The best example in our sample of an old (>10 Gyr), 
genuine globular cluster is the very metal-poor cluster 
NGC 2210 in the outskirts of the LMC (Fig. 5). It shows a 
richly populated RGB sequence, terminating abruptly at 
Mks = —5.95 ± 0.05 mag. There is no sign whatsoever of 
an AGB sequence, nor does it contain any IR objects. 
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Table 4. Ncar-IR photometry of cluster IR objects from 
the DENIS catalogue of the Magellanic Clouds (Cioni et 
al. 2000). 
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Table 5. New mid-IR photometry of cluster objects, ob- 
tained at a wavelength of 10.6 or 11.9 /an with TIMMI2 
at the ESO 3.6m telescope: flux densities are in Jy; the 
bands are listed by their central wavelengths in ^m. In 
addition, photometry in the MSX bands C, D and E are 
given where available. 



Cluster 



IR# 



10.6 



11.9 



12.1 14.' 



21.3 



NGC 1984 

BSDL1837 
NGC 1994 
HS 327-E 
KMHK 1603 



6.8±0.7 
0.39±0.08 

1.59±0.21 



5.28 9.9 19.8 



0.37: 
1.04 
0.44 



0.9 



1.1 



0.54±0.06 



The populous intermediate-age (~0.1-10 Gyr) clusters 
in our sample show an AGB sequence which generally con- 
tinues up to Mks ~ —8 mag and (J — K s ) ~ 1.8 mag (Fig. 
6). As mentioned earlier, the AGB sequences in the SMC 
clusters do not appear to be much bluer than those in the 
LMC. A good example is the SMC cluster NGC 419 of 
which the AGB is well populated up to (J — K s ) = 1.9 
mag, which is redder than some LMC clusters. NGC 1978 
is a good example where the RGB can be distinguished 
from the AGB: the stellar density is high up to Mk s = —6 
mag but then drops sharply as the RGB terminates and 
only the AGB continues: the AGB evolution is much more 
rapid than the RGB phase. Compared to NGC 2210 the 
RGB of more metal-rich intermediate-age clusters is red- 
der in J-K s . 

Many of the intermediate-age clusters contain one 
or more IR carbon stars (Fig. 6). The most extremely 
dust-enshrouded IR objects are much fainter in the K s - 
band than optically bright tip-AGB stars, as their dust 
envelopes become optically thick even at 2 /im (e.g., 
NGC419-IR2 at M Ks = -3.4 mag). Hot objects such as 
post-AGB objects and PNe are also fainter than tip-AGB 
stars at near-IR wavelengths, because they shine mostly 
at optical and UV wavelenghts. 



Table 6. Mid-IR photometry of cluster objects from 
ISOCAM: flux densities are in mJy; bands are listed by 
the central wavelength in /im. References are: l=Tanabe 
et al. (2004); 2=Trams et al. (1999b); 3=van Loon ct al. 
(2001a); 4=this work. 



Cluster IR# 


4.5 


6.7 


12.0 


14.3 


Reference 


Small Magellanic 


Cloud 










NGC 419 1 


71 


49 


37 




1 


2 


29 


56 


103 




1 


Large Magellanic 


Cloud 










HS33 1 






269 




2 


NGC 1783 1 


95 


75 


55 




1 


BSDL 1837 1 


212 




236 




4 


NGC 1978 1 


59 


50 


36 




1 


2 


61 


71 


71 




1 


4 


40 




15 




1 


HS 327-E 1 


187 


544 


433 


698 


3 


2 


15 


9.5 


5.7 


6.6 


3 


SL519 1 


118 


138 


155 


135 


4 


2 


1.0 


4.6 


2.6 


1.9 


4 


NGC 2121 1 


51 


34 


17 




1 



4.5 um| 



6.7 (jm | 



4 



S2 



IR2 



IR1 

S3 



SI 



N 




14.3 urn 



Fig. 3. Mid-IR images of SL 519 at 4.5, 6.7, 12.0 and 14.3 
fim, obtained with ISOCAM. 

The 25 Myr young SMC cluster NGC 330 displays a 
bright RSG branch at (J—K s ) = 1 mag, and massive Main 
Sequence stars at (J—K s ) — mag (Fig. 7). The older, 130 
Myr LMC cluster NGC 1866 also contains at least three 
very bright red giants. These must be massive (Mi > 4 
Mq) AGB stars. Such oxygen-rich M giants have near-IR 
colours that resemble those of red supergiants rather than 
the (redder) AGB carbon stars. The circumstellar dust 
shells of very luminous AGB stars and supergiants become 
optically thick only at much higher mass-loss rates than 
envelopes of lower luminosity AGB stars (cf. van Loon 
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et al. 1997, 1998). The dearth of reddened stars in young 
clusters does not however mean that no mass-losing giants 
are present. For instance, KMHK 292-IR1 and IR2 clearly 
show emission from circumstellar dust despite not being 
much reddened by it. 

NGC 1903 displays two sequences, separated in J-K s 
by a few tenths of a magnitude (Fig. 8). The bluer of these 
sequences probably traces a young population of massive 
AGB stars or red supergiants, whilst the redder sequence 
traces an intermediate-age AGB. The brightest star, IR3 is 
located near the cluster centre and seems to be associated 
with the bluer sequence. IR1 and IR2 are both located at 
the fringes of the cluster and, being carbon stars (van Loon 
et al., in preparation), must be associated with the older 
AGB sequence. Hence it is likely that NGC 1903 is young 
and that IR1 and IR2 are associated with an intermediate- 
age field population or a super-imposed cluster. This is 
confirmed in the spatial distribution of the stars (Fig. 9): 
the blue Main Sequence stars with (J — K s ) ~ mag and 
the "blue", luminous red giants with (J — K s ) < 1 mag 
or Mks < —8 mag form the cluster NGC 1903, whereas 
the redder, fainter red giants with (J — K s ) > 1 mag and 
Mk s > — 8 mag are concentrated towards the SW of the 
field. 

In NGC 1984 we see a similar superposition of a young 
cluster and an intermediate-age field population (Fig. 8). 
The bright red supergiant WOHG347 (IR2) and the 
highly obscured object IR1 are both located near the clus- 
ter centre. But IR3, located further out, is thought to be 
a PN and thus belongs to an intermediate-age population 
unrelated to the cluster. 



4.3. Cluster age determinations 

For the LMC clusters HS33, KMHK 292, KMHK 285, 
HS 270 and BSDL 1837 no age estimates could be found 
in the literature. We can make crude estimates of their 
ages on the basis of the IR colour-magnitude diagrams, 
where we adopt [Fe/H]= -0.4 (Table 1). BSDL 1837 is 
thus found to be of intermediate age. Because BSDL 1837- 
IR1 is with L ~ 10 5 L Q (see Table 8) twice as bright as 
an intermediate-age AGB star would ever be, it must be 
much younger than the cluster and therefore it is probably 
not a cluster member. 

For HS 33 we had already obtained Johnson V and 
Cousins I-band images with the 0.9m Dutch telescope 
at ESO La Silla, Chile, on New Year's Eve 1996. Multi- 
object photometry was obtained using DAOphot within 
MIDAS, and calibrated against standard stars observed 
in the SA98 and TPHE fields. The I versus V-I colour- 
magnitude diagram (Fig. 10) is compared with isochrones 
from Bertelli et al. (1994) for ages of 0.1, 0.2, 1 and 10 Gyr 
and a metallicity of [Fc/H]= —0.4. The cluster appears to 
have an age of t ^1301^0 Myr, populating much more of 
the upper Main Sequence around (V — I) ^0 mag than of 
the red clump around (V — I) <~1.0 and Mi <~ —0.3 mag. 
The >1 Gyr old RGB and AGB above the red clump are 



M 



1 1 1 1 i 1 1 1 1 i 

NGC 2210 



■I 



2 3 4 

J K s 



V 



2 4 
K R -L' 



Fig. 5. Example of an old cluster, NGC 2210, displaying 
an RGB terminating abruptly at Mk s = —6 mag (Section 
4.2). 
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Fig. 6. Examples of two populous intermediate-age clus- 
ters in the SMC and LMC, each with several dusty stars 
(Section 4.2). 



almost entirely devoid of stars in the direction of HS33, 
whilst these sequences and the red clump are heavily pop- 
ulated with field stars surrounding the cluster. HS 33- IRI 
appears to lie on the red extension of the 1 Gyr isochrone, 
but this is misleading as the object suffers from severe cir- 
cumstellar extinction and its colour and magnitude can be 
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Fig. 4. Cumulative colour-magnitude diagrams of absolute K s -band magnitude versus J-K s (top) and K s -L' (bottom), 
for the stars detected in the L' band in clusters in the SMC (open symbols) and LMC (filled symbols). 



traced back to the upper AGB of at most a few 100 Myr 
old. This is supporting evidence, both for the age of HS 33 
and for IRl's association with the cluster. 

4.4. Photometric classification of cluster IR objects 

The locations of dusty stars in IR colour-magnitude and 
colour-colour diagrams may clarify their nature. For in- 
stance, dust-enshrouded stars delineate a sequence in a 
J-H versus H-K s or K S ~L' versus J-K s diagram (Fig. 11). 
Oxygen-rich star are found at the blue end of this sequence 
and carbon stars at the red end (cf. Frogcl & Cohen 1982). 
The confirmed PNe, NGC 1852-IR1 and NGC 1984-IR3 lie 
in a disparate region of the near-IR colour-colour diagram, 
"below" the sequence in the J-H versus H-K s diagram (cf. 
Fig. 7 in Ferraro et al. 1995) and "above" the sequence in 
the K s -L' versus J-K s diagram. Whitelock (1985) shows 
that PNe can assume negative J-H colours if the He 1 1.083 



fim line is included in the J-band filter. This is the case 
in the SAAO filter suite but the 2MASS and ESO (SOFI) 
J-band filters do not include this line. 

In the K s - [12] and K s -[8.3] versus J-K s and versus 
K s -L' diagrams (Fig. 12), carbon stars and oxygen- rich 
red giants follow unique sequences chiefly as a function of 
optical depth. These sequences intersect each other and 
identification is therefore not always conclusive. We draw 
the sequences for oxygen-rich AGB stars (solid line) and 
carbon stars (dotted line) from van Loon et al. (1997) 
and Trams et al. (1999b), after small adjustments to ac- 
count for differences in photometric systems. Post-AGB 
objects stand out prominently in the K s -[12] versus J-K s 
diagram (van Loon et al. 1997) as they have strong dust 
emission but little extinction. Carbon stars also tend to 
occupy a distinct region in the [12]-[25] versus J-K s and 
versus K s -L' diagrams (Fig. 12) because their circumstel- 
lar envelopes are warmer and become optically thick more 
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Fig. 7. Examples of two populous young clusters; the 
SMC cluster NGC 330 and the LMC cluster NGC 1866 
(Section 4.2). 



easily, and when oxygen-rich envelopes become optically 
thick at mid-IR wavelengths the silicate absorption causes 
a reddening of the [12]-[25] colour. 

The IR diagrams suggest that NGC 1903-IR3, 
NGC1994-IR1 and NGC2100-IR1 are oxygen-rich stars. 
HS270-IR1, NGC1984-IR1, NGC1978-IR3 and SL482- 
IR1 have IR colours similar to the known low-excitation 
PNc NGC 1852-IR1 and NGC 1984-IR3, and might have 
a detached dust shell. 



5. Luminosities and mass-loss rates 

5.1. Modelling the spectral energy distributions 

The spectral energy distributions of the IR objects were 
modelled with the radiative transfer code dusty (Ivezic, 
Nenkova & Elitzur 1999). The density distribution is based 
upon a hydrodynamic computation of a dust-driven wind 
at constant mass-loss rate. The model was then scaled to 
match the overall observed SED, knowing the distances to 
the SMC and LMC objects, which then yields an accurate 
measurement of the bolometric luminosity. 

A blackbody was used to represent the underlying stel- 
lar radiation field, with a temperature of T c s = 2500 K in 
the case of carbon stars and T e g = 3000 K in the case of 
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Fig. 8. Examples of two LMC clusters, displaying a pecu- 
liar post-Main Sequence morphology (Section 4.2). 
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Fig. 9. Spatial distribution of L'-band detected sources 
near NGC 1903. The blue Main Sequence stars and "blue" , 
luminous red giants constitute the cluster NGC 1903, 
whereas the fainter, redder red giants are concentrated 
towards the SW. 



oxygen-rich, M-type stars unless the actual temperatures 
were known or the fit required a hotter star. 

For almost all objects the chemistry is known, and 
by default we chose amorphous carbon dust (Henning & 
Mutschke 1997) or astronomical silicate (Draine & Lee 
1984) for carbon-rich and oxygen-rich objects, respec- 




Fig. 10. Optical colour-magnitude diagram for the 
cluster (dots) and field (circles) around HS33 and 
IRAS 04496-6958 (IR1), with overlaid isochroncs for ages 
of 0.1, 0.2, 1 and 10 Gyr and metallicity of [Fe/H]= -0.4 
(Bertelli et al. 1994). 



tively. The extreme OH/IR objects NGC 1984-IR1 and 
HS 327-E-IR1 were better fit with warm, oxygen-deficient 
silicate from Ossenkopf, Henning & Mathis (1992), and the 
peculiar carbon-rich object NGC 1978-IR3 required the in- 
clusion of 20 per cent silicon carbide (SiC, from Pcgourie 
1988). A fit to SL 482-IR1 was obtained with cold, oxygen- 
rich silicate from Ossenkopf et al. (1992) and a very thin 
shell geometry to mimic an isothermal dust envelope at 
Idust = 750 K. We used either a single dust grain size, or 
a standard MRN grain size distribution (a power-law with 
exponent —3.5; Mathis, Rumpl & Nordsieck 1977). 

To obtain the mass-loss rate one has to know the dust 
grain density, for which we adopt p gra in = 3 g cm -3 , and 
the gas-to-dust mass ratio. Very little is known about the 
latter, which may vary according to the environmental 
conditions — metallicity and chemistry, radiation field, 
and evolutionary state. We adopt a reference value of 
tpQ = Pgas/pdust = 200 for solar metallicity, and scale 
this value according to the cluster metallicity (Table 1) as 

^ = ^ Q io-[ Fc / H i. 

The input parameters and output results are sum- 
marised in Table 8. Here we also list the initial mass for 
the cluster member descended from the heaviest Main- 
Sequence progenitor, M;, which is taken from the Bertelli 
et al. (1994) isochrones. These masses are representative 
for the cluster IR objects that are expected to be the most 
evolved stars in the clusters. For comparison we also list 
the Main-Sequence Turn-Off masses, defined here as the 
mass of the Main Sequence star that has reached its (first) 
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Fig. 11. Ncar-IR colour-colour diagrams to probe the na- 
ture of the stellar photospheres of the cluster IR objects. 

maximum in effective temperature as its core becomes de- 
pleted in hydrogen. 

5.2. Luminosities and nature of cluster IR objects 

Not surprisingly, all cluster IR objects are luminous (Fig. 
13b). Where the cluster has a RGB, the cluster IR objects 
were always more luminous than the RGB tip. The only 
exception is the faint object SL519-IR2, of which the ex- 
act nature is uncertain. Most objects straddle along the 
maximum luminosity reached within a cluster (the curve 
in Fig. 13b, derived from Bertelli et al. 1994). 

It is especially interesting to note that the luminous 
carbon star HS 33-IR1 matches very well the expected lu- 
minosity at the tip of the AGB for the cluster HS 33. This 
lends support to the interpretation of it having evolved 
from a massive AGB star progenitor (M; > 4 M Q ). 

All dust-enshrouded stars with progenitors of 1.3 < 
Mi < 2.2 M are carbon stars, whilst stars with progenitor 
masses of M\ > 6 M are always oxygen-rich. 

Our sample includes four clusters with Mi ~ 4 M Q . No 
IR-excess objects were found in the SMC cluster NGC 458, 
despite its richness and rather high metallicity. The LMC 
cluster objects SL453-IR1, HS 327-E-IR2 and SL519-IR1 
are carbon stars. These are the only cluster IR objects 



Table 7. List of cluster IR objects and their classifica- 
tion. Objects for which an L-band spectrum is available 
are referred to under "Lspec", where l=van Loon et al. 
(1999a); 2=Matsuura et al. (2002); 3=van Loon et al. 
(2003); 4=van Loon et al. (in preparation). Objects which 
are targets for the Spitzer Space Telescope are listed un- 
der "SST" by their SST programme number, where the 
Principal Investigators are J. Houck (#103), F. Kemper 
(#1094) and P. Wood (#3505). 



Cluster IR# Classification Lspec SST 

Small Magellanic Cloud 
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AGB carbon 
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AGB carbon 
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AGB carbon star 
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that have luminosities that are significantly lower than 
the maximum possible luminosity. However, the brightest 
IR object in these clusters, HS 327-E-IR1 is an oxygen- 
rich object of which the luminosity is right at the tip of 
the AGB for that cluster age. It has been argued (van 
Loon et al. 2001a) that the coincidence of carbon stars 
and more luminous oxygen-rich stars in populations with 
an age of t ~ 2 x 10 8 yr suggests that M; = 4 M Q is 
the threshold of Hot Bottom Burning (HBB; Boothroyd & 
Sackmann 1992), with the luminous oxygen-rich AGB star 
descendent from a progenitor with a slightly larger mass 
than this threshold and the carbon star having a slightly 
lower mass. Alternatively, a second parameter might be 
involved, such as rotation or a metallicity spread. 
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The late stages in the evolution of stars with progenitor 
masses in the range M\ = 5-8 M Q (shaded regime in Fig. 
13) are uncertain, with the Bertelli et al. (1994) models 
suggesting that the maximum progenitor mass of an AGB 
star is Mi ~ 5 M Q . The oxygen- rich red giant NGC 1903- 
IR3 is the only cluster IR object found in the three clusters 
in this mass range. It is more luminous than expected, 
which may be due to HBB (the Bertelli tracks do not 
include HBB). 

Unfortunately we did not find cluster IR objects in our 
sample around Mi ~ 10 M Q . Such stars are believed to 
become "super- AGB" stars that follow a thermal-pulsing 
AGB evolution but go on to ignite core-carbon burning 
(Ritossa, Garcia-Berro & Iben 1996). The supergiants of 
Mi = 14-19 M Q that we did find have luminosities that 
match very well the expected luminosities for the clusters 
they are associated with. 

5.3. The contribution of dust-enshrouded red giants to 
the integrated cluster IR brightness 

Dust-enshrouded giants do not contribute to the optical 
brightness of the cluster they are in, but they dominate 
the integrated cluster brightness at mid-IR wavelengths 
of A<;8 fim even in populous clusters. Thus, the mid- 
IR brightness of a cluster cannot be reliably predicted 
as the rare occurrence of a dust-enshrouded red giant 
is a stochastic event. However, if the combined effect of 
many clusters is considered, then a sufficient number of 
dust-enshrouded red giants will always be present and the 
fluctuations will be small. Computations of Single Stellar 
Populations (Brcssan, Granato & Silva 1998; Piovan, 
Tantalo & Chiosi 2003) may therefore yield reliable pre- 
dictions for elliptical galaxies but not for star clusters. 

With a central wavelength of 3.5 /zm, IRAC channel 1 
on-board the Spitzer Space Telescope is similar to the L- 
band; what is the contribution of the dust-enshrouded red 
giants to the cluster brightness in the L'-band? Consider 
the example of NGC 419, a populous cluster with two 
very red dust-enshrouded carbon stars. The combined 
brightness of these two IR objects is L' m = 8.66 mag. 
The combined brightness of all other stars in the ISAAC 
fields above the RGB tip [M^i = —6.4 mag; see Section 
5.5) is L' >RGB = 7.58 mag. Including all fainter stars 
down to i[ imit ~ 14 mag the cluster brightness becomes 
-^cluster = ^.49 ma g- This estimate excludes the (small) 
contribution of even fainter cluster members as well as 
stars outside of the ISAAC field of view. Thus, the two IR 
objects contribute < 30 per cent to the integrated cluster 
brightness. For smaller clusters the contribution can be 
much larger and the cluster L'-band brightness becomes 
completely stochastic. 

5.4. Mass-loss rates of cluster IR objects 

Most cluster IR objects are fairly evenly distributed over 
more than a decade in mass-loss rate with a maximum of 



M ~ 10~ 4 M Q yr _1 (Fig. 13c). The spread in mass-loss 
rates for a given progenitor mass is due to evolutionary 
effects, as the stars are not always captured at the highest 
mass-loss rate. This effect can be traced by the position 
of the star in the Hertzsprung-Russell diagram, and van 
Loon et al. (1999b, 2005) indeed show empirical evidence 
for a dependence of the mass-loss rate on the bolomctric 
luminosity and stellar effective temperature. 

However, upon closer inspection it seems that the lower 
mass progenitors, Mi ~ 1.4 M , reach mass- loss rates that 
are lower than those of somewhat higher mass progenitors, 
Mj ~ 2 M , whilst the highest mass-loss rate occurs at 
an even higher progenitor mass of Mi ~ 4 M Q . A similar 
tendency is seen amongst the massive progenitors of Mi > 
13 Mq. As both the (maximum) luminosity and cluster 
metallicity increase with increasing progenitor mass, it is 
difficult to ascertain the physical cause for the dependence 
of the mass-loss rate on progenitor mass. 

Objects with detached shells (triangles in Fig. 13) are 
always the stars with the highest mass-loss rate at a given 
progenitor mass. For these stars the mass-loss rate refers 
to the detached dust shell and is thus a historic account of 
the mass loss experienced some time in the past (typically 
a few centuries ago). Hence the most likely explanation 
is that these are stars which have just emerged from the 
most intense superwind phase in their evolution. Amongst 
these, the luminous OH/IR object NGC 1984-IR1 had an 
extremely high mass-loss rate of M ~ 4 x 10~ 3 M Q yr _1 . 
Modelling the SED using a cool star of 2500 K hardly 
affects the estimate of its luminosity but the mass-loss 
rate would be even higher by ^50 per cent. 

No cluster IR objects are found in clusters with cither 
M < 1.3 M or [Fe/H]< -0.9 (Fig. 13a). The L'-band 
observations were sensitive enough to have detected all 
AGB stars and stars on the top magnitude of the RGB 
even if they had no L'-band excess, and the old clusters 
included some populous globular clusters with plenty of 
red giants. The threshold for identification as an IR object 
is M m ; n ~ 10~ 6 M Q yr _1 (Table 8). Lower mass-loss rates, 
a lower dust content or a shorter duration of the mass- 
loss episode (or a combination of these) might explain the 
dearth of cluster IR objects with M ; < 1.3 M Q . 

Unfortunately, there is an age-metallicity degeneracy 
in the sense that there are no old metal-rich clusters or 
young metal-poor clusters. Galactic globular clusters, that 
are all old but span a metallicity range from [Fc/H]< —2 to 
~solar, show a clear threshold around [Fe/H]^ —1 below 
which no large-amplitude (Mira-type) pulsating stars are 
found (Frogel & Elias 1988). With pulsation believed to be 
an essential first stage in the mass-loss of red giants, failure 
to pulsate strongly enough might explain the absence of 
dust-enshrouded red giants at [Fe/H]< —0.9. 

5.5. Timescales of mass-loss 

To estimate the impact of mass loss in the superwind stage 
relative to the total stellar mass, we could compare the 



van Loon, Marshall, Zijlstra: Dust-enshrouded giants in clusters in the Magellanic Clouds 



15 



Table 8. Initial metallicity, [Fe/H], and mass, Mi, and Main-Sequence Turn-Off mass, Mto, and the input and results 
of the modelling with dusty: stellar effective temperature, T e fj, temperature at the inner radius of the dust envelope, 
Tdustj dust type (Sil-DL=astronomical silicate, Draine & Lee 1984; Sil-Oc=cold oxygen-rich silicate, Ossenkopf et al. 
1992; Sil-Ow=warm oxygen-deficient silicate, Ossenkopf et al. 1992; AmC=amorphous carbon, Henning & Mutschke 
1997; SiC=silicon carbide, Pegourie 1988), grain size, a (where a range is given, these represent the minimum and 
maximum grain sizes for a standard MRN distribution — Mathis et al. 1977), bolometric luminosity, L, and total 
(gas+dust) mass-loss rate, M. 



Cluster IR# 


[Fe/H] 


M 


M TO 


T e s 


Tdust 


dust 


a 


log(L) 


log(M) 


cluster 








(M Q ) 


(M ) 


(K) 


(K) 


type 


(H 


(Le) 


(M /yr) 


member 


Small Magell 


anic 


Cloud 




















NGC419 


1 


-0.60 


1.98 


1.57 


2800 


1000 


AmC 


0.2 


4.00 


-5.00 


yes 




2 








2800 


700 


AmC 


0.1 


3.85 


-4.34 


yes 


Large Magelh 


%nic 


Cloud 




















HS33 


1 


-0.4 


4.8 


4.0 


2500 


720 


AmC 


0.01-0.1 


4.61 


-4.54 


yes? 


KMHK292 


1 


-0.4 


15.0 


11.9 


3000 


400 


Sil-DL 


0.1 


5.37 


-4.62 


yes 




2 








3000 


350 


Sil-DL 


0.1 


5.17 


-5.20 


yes 


KMHK285 


1 


-0.4 


2.21 


1.73 


2800 


600 


AmC 


0.01-0.1 


3.99 


-4.57 


yes 


NGC 1783 


1 


-0.75 


1.32 


1.16 


2800 


1100 


AmC 


0.1 


3.93 


-4.90 


yes 


NGC 1852 


1 


-0.85 


1.31 


1.16 


28600 


500 


AmC 


10 


4.29 


-4.25 


yes 


NGC 1903 


1 


-0.4 


6.3 


5.3 


2800 


1300 


AmC 


0.01-0.1 


3.99 


-5.11 


no 




2 








2800 


1500 


AmC 


0.01-0.1 


3.84 


-5.43 


no 




3 








3000 


1200 


Sil-DL 


0.01-0.1 


4.66 


-4.87 


yes 


HS270 


1 


-0.4 


2.21 


1.73 


2500 


225 


AmC 


0.01-0.3 


3.94 


-4.15 


yes 


SL453 


1 


-0.4 


4.0 


3.3 


2500 


1100 


AmC 


0.01-0.1 


3.70 


-4.75 


yes 


SL482 


1 


-0.4 


13.7 


11.1 


20000 


750 


Sil-Oc 


0.1 


5.19 


-4.02 


yes 


NGC 1984 


1 


-0.90 


19.1 


14.3 


15000 


400 


Sil-Ow 


0.01-0.1 


5.39 


-2.38 


yes 




2 








4000 


1000 


Sil-DL 


0.1 


5.10 


-4.25 


yes 




3 








31000 


650 


AmC 


0.1-0.2 


4.47 


-4.91 


no 


BSDL1837 


1 


-0.4 


2.21 


1.73 


2500 


1000 


Sil-DL 


0.1 


5.04 


-5.03 


no 


NGC 1994 


1 


-0.24 


14.1 


11.7 


2500 


600 


Sil-DL 


0.1 


4.90 


-4.02 


yes 


NGC 1978 


1 


-0.66 


1.42 


1.30 


2500 


1000 


AmC 


0.1 


3.73 


-5.05 


yes 




2 








2500 


1100 


AmC 


0.1 


3.77 


-4.80 


yes 




3 








2500 


300 


AmC+SiC 


0.01-5.0 


3.70 


-4.39 


yes? 




4 








2500 


1000 


AmC 


0.1 


4.19 


-5.46 


yes 


HS327-E 


1 


-0.4 


4.0 


3.3 


2500 


800 


Sil-Ow 


0.01-0.1 


4.45 


-3.90 


yes 




2 








2500 


200 


AmC 


0.1 


3.89 


-5.83 


yes 


SL519 


1 


-0.4 


4.0 


3.3 


2500 


650 


AmC 


0.01-0.1 


3.84 


-4.67 


yes 




2 








2500 


830 


AmC 


0.2 


2.58 


-6.26 


yes? 


NGC 2100 


1 


-0.32 


14.3 


11.7 


2500 


300 


Sil-DL 


0.1 


4.87 


-4.69 


yes 


NGC 2121 


1 


-0.61 


1.31 


1.15 


2500 


1200 


AmC 


0.2 


3.95 


-5.36 


yes 


KMHK 1603 


1 


-0.60 


2.17 


1.68 


2500 


600 


AmC 


0.1 


4.09 


-4.33 


yes 



Table 9. Total cluster mass, mass-loss rate and feedback timescale, and number of cluster IR objects, number of stars 
brighter in L' than the RGB tip, their ratio and the cluster mass per super-RGB-tip or IR star, for different initial 
mass or cluster age. 
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total mass-loss rate with the total cluster mass (Table 
9). We do this for different ranges in the initial mass of 
the mass-losing stars, corresponding to different ranges in 
cluster ages (these are approximate as they also depend 
somewhat on initial metallicity) . The older clusters in our 



sample comprise an order of magnitude more mass than 
the younger clusters, but the total mass- loss rate is much 
higher in the younger clusters. This results in a very strong 
age dependence of the timescale for the cluster mass to de- 
crease, t mass loss = Y2 M clu stor/ Y2 M - 
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Fig. 13. Luminosities and mass-loss rates of cluster IR objects as a function of their initial mass. Panel (a) displays 
all clusters in a mass-metallicity diagram. The mass corresponds to the initial mass of the currently most evolved star 
in the cluster. Oxygen-rich giants, carbon stars and objects with detached shells are represented by discs, boxes and 
triangles, respectively, where open symbols indicate that the IR object is most likely not a physical member of the 
cluster. The curve in the luminosity diagram (b) traces the maximum luminosity achieved, from Bertelli et al. (1994) 
models of [Fe/H]=— 0.4 and —0.7 for masses larger and smaller than 2 M©, respectively. The shaded area indicates the 
regime between Mi = 5 and 8 M Q where the transition between AGB and RSG is expected. In panel (c) the dashed 
horizontal at M = 10~ 4 M Q yr _1 is for reference only. 



In principle this timescale is the time it takes for the 
entire cluster mass to be dissipated if the current mass-loss 
rate were sustained. For stellar systems in general (such 
as galaxies) it is the timescale for the feedback mechanism 
to return all mass back into the ISM. It is remarkable that 
although the estimated timescales span a range of more 
than two orders of magnitude, they are always about 2 to 
3 times the cluster age. 

Another approach to quantify the timescale for mass 
loss is to compare the number of cluster IR stars to the 
number of stars in an evolutionary stage for which the 
timescale is known. For the latter we choose to count the 



number of stars in our L'-band images that are brighter 
than the RGB tip, which occurs at Ml> = —6.4 mag 
(Table 1). This has a different meaning for clusters of dif- 
ferent ages, for instance young clusters do not contain any 
RGB stars, and will need to be calibrated against stellar 
evolutionary models. But for a wide range of cluster age 
it does provide a rough measure for the number of AGB 
stars on the thermal-pulsing AGB, whilst for the younger 
clusters it provides a measure for the number of post-Main 
Sequence supergiants. The latter is a consequence of our 
definition of the RGB tip in the L'-band which, coupled 
with large bolometric corrections for blue Main Sequence 



van Loon, Marshall, Zijlstra: Dust-enshrouded giants in clusters in the Magellanic Clouds 



17 



stars discriminates between massive Main Sequence stars 
and redder post-Main Sequence supergiants. 

The ratio of IR objects to the total number of stars 
above the RGB tip is the same within a factor two for 
all clusters except the oldest clusters, in which we did not 
find any IR object and which anyhow contain very few 
stars above the RGB tip. The ratios indicate that the su- 
pcrwind timescales are 10-20 per cent of the time stars 
spend at Ml* < —6.4 mag. There is a hint of a trend 
for this timescale to be longer for older clusters (up to 20 
per cent) compared to younger clusters with AGB stars 
(15 per cent) and the youngest clusters with red super- 
giants (10 per cent). Intermediate-mass stars spend about 
10 6 yr above the RGB tip, which implies a timescale of 
1-2x10 yr for the superwind phase in which the stars 
become bright IR objects. This corresponds to only a few 
thermal pulses, presumably the final ones, and is consis- 
tent with the fact that the IR objects are found within a 
magnitude from the AGB tip: this magnitude range can 
be fully accounted for by luminosity variations over the 
thermal pulse cycle of a factor ^thrce (the evolutionary 
timescale for an AGB star to increase in luminosity by a 
magnitude is about 10 6 yr). These timescales and timing 
at the end of AGB evolution confirm model predictions 
(e.g., Vassiliadis & Wood 1993). 

For the low- and intermediate-mass AGB stars one can 
estimate the importance of the superwind phase compared 
to the more moderate AGB mass loss. Stars of initial 
masses 1.3-1.8 M will loose about 1 M Q on the AGB. 
We find 36 stars above the RGB tip in the cluster sample 
(Table 9), which should therefore loose ~36 M Q in total. 
The ratio of this mass lost and the integrated mass-loss 
rate of 1.4 x 10~ 4 M Q yr -1 yields a duration of 3 x 10 5 
yr required for this mass to be lost exclusively during the 
superwind phase. The same exercise for stars in the 1.8-3 
M Q range, which loose about 2 M Q on average, yields an 
identical estimate for this timescale. Given the above de- 
rived superwind timescale of 1-2 xlO 5 yr this implies that 
a significant fraction of the mass loss (30-70 per cent) in- 
deed occurs during the superwind phase, again confirming 
model predictions (e.g., Wachter et al. 2002). 

The underlying cluster mass associated with each 
super-RGB-tip star increases as the cluster ages. This can 
be turned around and used to estimate the underlying 
cluster mass when the brightest (in L') cluster members 
can be counted, provided that the cluster age is known. 
We thus obtain a first crude estimate of such relation: 

log y Mg)] = 3.6(±0.2) + 0.30(±0.06) i[Gyr]. (3) 

N >RGB 

6. Summary of conclusions 

We present the results of a search for dusty evolved stars 
in star clusters in the Magellanic Clouds. The survey was 
performed in first instance through imaging in the near- 
IR V band (at 3.8 /an) with ISAAC at the ESO-VLT. 
Comparison with imaging photometry at shorter (J, H 



and/or K s ) and longer (MSX, IRAS, ISO and ground- 
based mid-IR) wavelengths allowed the identification of 
objects that are surrounded by dust causing extinction in 
the optical and near-IR and excess emission in the ther- 
mal IR. Targets were selected to include populous clusters 
spanning a range in age and metallicity, and clusters that 
were already known or suspected to be associated with 
sources of bright IR emission. 

Out of 9 clusters in the SMC and 29 clusters in the 
LMC, 19 clusters (of which one in the SMC) were found 
to contain a total of 30 stars with IR excess emission. 
Of these, 4 stars are probably not physically associated 
with the cluster. We establish the nature of the major- 
ity of the cluster IR objects (i.e. their spectral class and 
evolutionary state). They are all highly evolved, and in- 
clude post-superwind objects such as Planetary Nebulae. 
No dusty stars are found for initial masses M\ < 1.3 M Q or 
initial mctallicities [Fe/H]< —0.9. All objects in the range 
1.3 < Mi < 2.2 M Q are carbon stars, and all objects with 
M; > 6 M Q are oxygen-rich giants or supergiants. 

We estimate the bolometric luminosities and mass-loss 
rates by modelling the spectral energy distributions with 
the DUSTY code. The IR objects are the most luminous 
cluster objects, nearly always at the maximum luminosity 
expected for their initial mass and metallicity. The mass- 
loss rate increases with larger progenitor mass, which may 
be due to a dependence on the initial metallicity or simply 
the luminosity. Post-superwind objects always have the 
highest mass-loss rates associated with them, where these 
mass-loss rates refer to the previous superwind stage. 

Stars brighter than the RGB tip in the L'-band spend 
10-20 per cent of their time (~ 10 5 yr) as IR objects with 
mass- loss rates that exceed a few 10~ 6 M Q yr -1 . About 
half of the mass lost by low- and intermediate-mass AGB 
stars occurs during this superwind phase. 

The number of stars above the RGB tip (in L') can be 
used to estimate the cluster mass, a relation for which we 
present a first crude estimate. 
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